Theory for supersolid 4 He 
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Although both vacancies and interstitial have relatively high activation energies in the normal 
solid, we propose that a lower energy bound state of a vacancy and an interstitial may facilitate 
vacancy condensation to give supersolidity in 4 He . We use a phenomenological two-band boson 
lattice model to demonstrate this new mechanism and discuss the possible relevance to the recently 
observed sup erfluid- like, non-classical rotational inertial experiments of Kim and Chan in solid 4 He. 
Some of our results may also be applicable to trapped bosons in optical lattices. 
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Recently Kim and Chan have reported observation of 
superfluid-like, non-classical rotational inertial (NCRI) 
behavior in solid 4 He, both when embedded in Vycor 
glass pj and in bulk 4 He 0)- Their experiments have re- 
vived great interest in supersolidity with both crystalline 
and superfluid orderings in helium. The possibility of 
a supersolid phase in 4 He was theoretically proposed by 
Chester 3], Leggett J], Saslow || and by Adreeve and 
Lifshitz [3 in 1970's. Adreeve and Lifshitz proposed Bose 
condensation of vacancies as the mechanism for superso- 
lidity. Chester speculated that supersolidity cannot exist 
without vacancies and/or interstitials 3], a claim made 
more rigorous recently by Prokofev and Svistunov 0. 
Experiments and more sophisticated microscopic calcu- 
lations have, however, provided constraints to any theory 
for supersolid in 4 He. The NMR experiments 0, 0] on 
solid 3 He rule out a non-negligible zero point vacancy 
concentration. The energy of a vacancy in solid 4 He was 
estimated to be about 10K by the x-ray scattering exper- 
iment [Rj and to be 15K in a theoretical calculation [Tl|. 
The energy of a pure interstitial state has recently been 
calculated to be about 48 ± 5K 12]. Assuming that the 
observed supersolidity is a genuine bulk phenomenon, the 
quandary is thus how defects with such high activation 
energies can condense at low temperature. 

In this Letter, we propose a possible solution to this 
quandary. In addition to vacancies and interstitials, there 
is a third type of defects, with relatively low excitation 
energy, which corresponds to a bound state of a vacancy 
and an interstitial, henceforth called an " exciton" . While 
such excitons do not carry mass and will not contribute 
to supersolid phenomena like those observed by Kim and 
Chan [HQ, they can facilitate supersolidity by two mech- 
anisms. First, vacancies can condense above a back- 
ground of excitons easier than above the defect free (DF) 
normal state, so that the condensation energy can more 
than compensate for the exciton excitation energy. Sec- 
ond, in a background formed of a coherent mixture of 
the DF state and the exciton, the effective kinetic en- 
ergy of vacancies and interstitials can be enhanced due 
to constructive interference between hopping processes 
involving the DF state and those involving the exciton. 



The essence of our theory is that while, consistent with 
all known experiments, the normal solid state is the DF 
state, the supersolid state results from condensation of 
vacancies and/or interstitials about a defect rich back- 
ground of excitons. This physics is shown quantitatively 
using a phenomenological two-band lattice boson model 
to represent the defects in solid 4 He. Using mean field 
theory (MFT), we show that superfluidity in solid he- 
lium can exist in parameter regimes qualitatively con- 
sistent with all the known experiments and microscopic 
calculations. We will argue that the key results will hold 
beyond MFT, and indeed are rendered more robust by in- 
clusion of quantum fluctuations. Because of the "vacuum 
switching" between the normal and supersolid states, the 
transition at zero temperature is generically first order. 

We start with the lattice as defined by the periodicity 
of crystalline helium. On each lattice site, we consider 
two single-particle localized states. The lower energy 
state (a-state), with energy — e a , has its maximum on 
the lattice site. The other state (6-state), with a higher 
energy — e a + A, is less localized and has maxima dis- 
tributed with hep symmetry away from the lattice site. 
Because of the strongly repulsive cores, each of this state 
can hold at most one He atom, i.e. for each state, helium 
behaves as hard core bosons. On the other hand, because 
of the spatial separation between the a and b states, an 
atom in the o-state repels one in the 6-state with a large 
but weaker strength U. An atom on one site can tun- 
nel to a neighboring site from a to a-state, b to b -state, 
and a to 6 state with hopping amplitudes t a , t^, and t a b 
respectively, taken all to be real and non-negative. In 
the DF state, we have one helium atom occupying the 
a-state on each lattice site. Due to the hard core condi- 
tion, the 4 He atoms in this state are immobile, and this 
state is a normal solid, which is consistent with the as- 
sumption that defects are necessary for supersolid. To 
study defects, it is convenient to consider the DF state 
as the defect vacuum state, from which vacancies and in- 
terstitials can be created. The term vacuum will always 
refer to the DF state henceforth in this article. Let a\ 
be the creation operator for the vacancy by removing a 
helium atom from this vacuum state at the site i, and b\ 
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the creation operator for the interstitial by creating an 
atom to the 6-state on the site. Note that the interstitial 
defined this way can be viewed as a quantum generaliza- 
tion of the classical interstitial. The hard core conditions 
specified above implies these operators are also hard core 
boson operators. Our model Hamiltonian for defects in 
solid helium is then, 

H = ^2 e a.nj,a + tbnj.b — Un^aTlj.b 
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In the above equation, 
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(previously defined) and are the energy cost to create 
a helium vacancy and an interstitial respectively. The 
repulsion U between He atoms translates into an attrac- 
tion of strength U between a vacancy and an interstitial 
on the same site. The energy of a localized exciton is 
A = e a + e b — U . For simplicity, in what follows we will 
consider t ab = 0. 

To illuminate the effects of excitons, let's first discuss 
the non-interacting case U = 0. In this limit, the va- 
cancies and interstitials are decoupled. For e b — > oo, 
our model is reduced to the vacancy model proposed by 
Andreeve and Lifshitz |g. Let z(= 12) be the number 
of the nearest neighbors in a hep lattice, the onset of 
Bose condensation of vacancies (interstitials) is given ex- 



actly by e 



a(b) 



zf. 



0(6) 



, which coincides with hav- 



ing zero activation energy for a single defect, a condition 
that is not supported by experiments ^3] an d theoretical 
estimates [ill ll^ for 4 He. This is the quandary we stated 
in the introduction. 

We now consider the case U > 0. In this case, a va- 
cancy and an interstitial tend to bind together to form a 
local exciton. If the exciton energy A is small (large U), 
the presence of the excitons in the ground state enhances 
the kinetic energy of the vacancy (or interstitial) , which 
may lead to the condensation of these defects. Since we 
are primarily interested in the large U case, it is impor- 
tant to treat the on-site attractive interaction accurately. 
In order to do so, we use the single-site mean field ap- 
proximation (MFA) by decoupling the kinetic terms as 
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where the spatially uniform Bose condensation order pa- 
rameters a =< a,i > and b =< bi > are determined 
self-consistently. The single site mean field Hamiltonian 
for the hard core bosons are then solved exactly. This 
MFA gives the correct exact conditions for onset of su- 
perfluidity at U = 0. 

The ground state phase diagrams obtained within the 
MFA in the parameter space t a and t b are shown in Fig. 
1(a) for e a = e b = 1, U = 1.9 and in Fig. 1(b) for e a = 1, 
e b = 4, and U = 4.8. They represent a more symmet- 
ric and a strongly asymmetric cases respectively. Within 
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FIG. 1: (a) Mean field theory phase diagram of Hamiltonian 
(1) in the symmetric case with e a = £b = 1 and U = 1.9. V-SF 
(I-SF): vacancy (interstitial) superfluid phase, (b) Same as in 
(a) in the asymmetric case e a = 1, e& = 4 and U = 4.8. See 
Fig. 2 for the snapshots of these phases in real space. (c)The 
variational energy of the symmetric model with e a = &b = 1 
and U = 1.9 as the function of the SF order parameter for 
t — t a — tt — 0.056 (normal solid phase), t — 0.06 (at the 
phase boundary) and t — 0.064 (VI-SF phase), (d) Same as 
in (c) with t a =tb = 0.06 for various values of U. 



the MFA, we found five different phases, characterized by 
the order parameters a and b together with the vacuum, 
vacancy, interstitial, and exciton defect densities no,ny, 
ni, and n ex . They are (1) The normal DF solid phase. 
(2) The vacancy superfluid phase (V-SF(A)), where only 
the vacancies condense (a ^ 0) above the DF background 
(nj = n ex = 0). This phase is the same as the vacancy 
state of Andreeve and Lifshitz. (3) The corresponding in- 
terstitial superfluid phase (I-SF). (4) An alternative va- 
cancy superfluid phase (V-SF(B)), where the vacancies 
condense above a background of excitons. (5) The VI- 
SF phase, where we have both vacancies and interstitials 
condensing above a background of a coherent mixture of 



3 



Normal Solid 



























-•- 

1 


2 


-•- 

3 


4 


5 


6 


7 


8 




V 


v 


12 


V-SF(A) 


























1 


• 

2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


12 


l-SF 






















• 




1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


12 



V-SF (B) 



























1 


-•- 

2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


12 



VI-SF 



FIG. 2: Schematic illustration of snapshots of five phases 
shown in fig^a) and (b) in terms of the occupation of he- 
lium atoms at the lattice sites. A lattice site is represented 
by a numbered block. A solid dot represents a helium atom. 
At each site, emptiness represents a vacancy, occupation at 
the lower (higher) level represents a vacuum (exciton), and 
occupation at the both levels represents an interstitial. 



the vacuum and the exciton. A snapshot of each of these 
phases in terms of the occupation of helium atoms in real 
space is illustrated in Fig. [21 

To appreciate the phase diagrams and the new physics 
arising from the interaction U, we show the phase bound- 
aries of the non-interacting model (U — 0) by dashed 
lines in Fig. 1(a) and (b). As remarked earlier, these 
boundaries coincide with the vanishing of the activation 
energy of an isolated vacancy or interstitial. Since U 
plays no role if only vacancies or only intersititials are 
present, these are also the phase boundaries for normal 
solid to V-SF(A) or FSF transition even for U > 0. Thus, 
for such SF states, their corresponding normal (i.e. un- 
condensed) states are not the DF state but contain a 
finite density of defects, which is not the case for super- 
solid 4 He. The most striking of our results is that the 
interplay between vacancies (and/or intersitials) and ex- 
citon defects can lead to SF even in parameter regimes 
where the normal state is stable against the generation of 
uncondensed defects (hatched areas in the figures). We 
propose this as the reconciliation between activated de- 
fect behavior at high temperature (T) and supersolidity 



at low T in 4 He. In the hatched regime, the DF state 
is metastable but not the global minimum energy state. 
In Fig. 1(c) and 1(d), we illustrate this by plotting the 
variational energy E(a) based on the MFA for the 
symmetric model with t a = = t, so that a = b, with 
increasing t (Fig. 1(c)) and U (Fig. 1(d)). The existence 
of two minima are clearly seen. With increasing t or U, 
the Bose condensed state becomes lower in energy than 
the normal state in the hatched areas. The transition is 
first order due to "vacuum switching": the normal state 
at a — is the DF state with no excitons (n ex = 0), 
while the background that vacancies and intersitials con- 
dense above to form the VI-SF is the one with a finite 
density of excitons (n ex ^ 0). In what follows we discuss 
in more details the physics of the asymmetric and sym- 
metric cases, focusing on the parameter regimes of the 
hatched areas. 

Due to the high He atom density and strong repul- 
sive cores, solid 4 He should fit the very asymmetric case 
with £ft/e a >> 1 in our model (Fig. 1(b)). Because 
the 6-state is less localized compared to the a-state and 
because of the exponential dependence of the overlap in- 
tegral on state size, we expect % » t a . For small t a , 
the phase in the hatched region is the V-SF(B) phase. 
To understand how this phase comes about, we consider 
first t a = 0. Here, the activation energies of a vacancy 
and an exciton are e a and A respectively independent of 
tf,, so the DF state is locally stable. Since t a — 0, vacan- 
cies cannot hop, and there is no possibility of vacancy 
condensation above the DF state. However, if we take 
a background of excitons on every site, then vacancies 
can now hop with amplitude and the exciton state 
can be unstable with respect to a vacancy condensation 
that gives 5^0. The onset of this V-SF(B) instability 
is ztb = e a — A. For % greater than this value, E(b) has 
a double minima behavior similar to that shown in Fig 
1(c) for the symmetric case, with E(b — 0) initially the 
lower energy. However, as t\> increases (but still less than 
the value necessary for spontaneous generation of inter- 
stitials), the vacancy condensation energy can become 
large enough to overcome the required vacancy and exci- 
ton activation energies given by nye a + (1 — ny)A, and 
a first order transition from the normal DF state to the 
V-SF(B) state occurs. Note that the transition is accom- 
panied by a "vacuum switching", in that the vacancies 
are condensing not above the DF state, but above the ex- 
citon state. The V-SF(B) has 5^0, and ny + n ex = 1. 
As t a increases, it becomes advantageous to mix in some 
vacuum component to allow vacancy hopping through 
t a , so that both a and b are non-zero. Depending on the 
value e&, some interstitial condensation will also occur. 
The V-SF(B) phase then makes a second order transi- 
tion into the VI-SF, characterized by both a and b ^ 0, 
and ny, n ex all ^ 0. We will discuss this phase further 
in the symmetric case, where it will feature more promi- 
nently. In this strongly asymmetric case, whether we 
take the SF state to be the V-SF(B) or the VI-SF state, 
we have ny » nj. Thus, a prediction of our theory is 
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that while the normal state is a commensurate solid, the 
supersolid will have incommensurate density, which can 
be confirmed by neutron scattering experiments in the 
supersolid phase. Furthermore, because of the presence 
of excitons in the supersolid (i.e. the He atom resides in 
the less localized b— state rather than the a— state, or 
some linear combination of the two), the local He den- 
sity in a unit cell will change with the transition into the 
supersolid. This can be confirmed from the form factor 
of neutron scattering or with a local probe. 

Next we look at the symmetric case. Although this 
probably does not describe solid 4 He, it may be applica- 
ble to a system of trapped bosons in an optical lattice, 
where the periodic potential is imposed externally rather 
than internally generated. In such a system, the ratio 
£b/e a can be tuned by tuning the optical potential. The 
interesting phase here is the VI-SF phase in the hatched 
region. Again, the key question is what causes the addi- 
tional non-trivial SF solution since neither vacancy nor 
the interstitial alone can condense at T = 0. Unlike the 
asymmetric case, where nj in this phase basically plays 
no role, and indeed — > as — > oo, here, no, ny, nj, n ex 
are all non-negligible in this SF phase. Within the MFA, 
the eigenstates are direct product of single-site states. 
In this ground state, the state on each site is a coher- 
ent mixture of the DF state, the vacancy, the intersti- 
tial, and the exciton. This coherence allows constructive 
interference between the various hopping processes, thus 
enhancing the effect of kinetic energy over that when only 
vacancies or interstitials are present or when they hop in 
a background of only DF states or only exciton states. 
This coherent effect is best understood by examining the 
highly symmetric case of = e a = e, ti, = t a = t and 
U = 2e a (so that A = + ). By symmetry, we also have 
a = b. Because of the symmetry between the vacancy 
and the interstitial, and between the vacuum and the ex- 
citon states, the MF energy Emf{o) can be analytically 

found to be | — J (|) 2 + (2zta) 2 + 2zta 2 , from which 
the condition for the SF state is found to be e — 2zt < 0. 
Compared with the SF condition for vacancy only con- 
densation above the DF background, we see that the ef- 
fective hopping integral increases from t to 2t due to the 
coherent mixture of all 4 states. For this artificial highly 
symmetric case, the normal to VI-SF transition is second 
order due to A being infinitesimal. As A increases, the 
transition becomes first order. This is what is shown in 
Fig lc as the transition is crossed by increasing t. Alter- 
natively, the transition can be crossed by increasing U, 
as shown in Fig 1(d). 

Although our results so far are obtained using the 
MFA, we believe the central conclusion, that supersolid 
can occur even when defects like vacancies and intersti- 
tials have relatively high activation energies is correct. In 
the lattice model of (1), the excitations in the supersolid 
phase are gapless and the excitations in the normal solid 
are gapful. Therefore, the quantum fluctuations are ex- 
pected to further stabilize the supersolid phase. In our 



MFT, there are four supersolid phases in terms of the 
order parameters a and b. However, there is no differ- 
ence between them in the type of off diagonal ordering 
of the underlying 4 He atoms. Thus, they are not really 
distinct phases in the sense of different symmetry states, 
but rather differ in the physics behind the condensation 
as discussed above. Indeed, if fluctuations are included 
or if we consider non-zero t a b, then the sharp distinction 
between these "phases" becomes crossover. 

Following Leggett's paper £|, we can derive the SF 
density detected by the NCRI experiment, which is 

Ps = Pam 0/ m* + p b m 0/ ml, (3) 

where m* and to£ are the effective mass of the vacancy 
and interstitial bands respectively determined by the 
hopping integral t a and i& and the lattice structure, and 
mo is the bare mass of 4 He. For the purpose of illus- 
tration, we will assume that to£ — mo at the value of 
<5 = 1/3. Because the transition at T = is first order, 
the SF density has an abrupt jump at the criticality. For 
reasonable values of parameters e a — 1, t\> = 4, t a = 0.07 
and ti, = 0.2, we find the jump to be ~ 9%, which is about 
an order larger than the SF density measured in Kim and 
Chan's experiments [2}. The discrepancy is partly due to 
the simplicity of the model. For instance, the critical SF 
density ps is found to be ~ 3% by including an on-site 
interband hopping t® b — 0.12. Furthermore, the quan- 
tum phase fluctuations mentioned above will decrease ps 
at criticality by the dual effects of directly decreasing the 
order parameters and by reduction of the critical value U c 
due to stabilizing of the SF phase relative to the DF state. 
We also note that the experiments of Kim and Chan are 
performed on granular rather than single crystal, so that 
the observed ps may be governed by Josephson effect and 
therefore can be considerably smaller than the intrinsic 
value of homogenous supersolid 4 He. 

We have seen that the T=0 transition should be first 
order. A careful examination of the finite T transition 
should include also the effect of phase fluctuations. The 
finite T transition as well as the collective excitations in 
the various SF phases of our model will be discussed in 
a future publication. 

In conclusion, we have proposed a solution to how 4 He 
can become a supersolid at low T when both vacancies 
and interstitials have relatively high activation energies 
in the normal solid. In our theory, the presence of low en- 
ergy bound vacancy-interstitial defects facilitate the con- 
densation of vacancies. In this theory, the normal state 
is a commensurate solid while the supersolid is incom- 
mensurate. Furthermore, the local helium density in the 
unit cell has different profiles in the two phases. These 
predictions can be tested experimentally. 
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our attention to this interesting topic as well as many 
stimulating discussions. 



5 



[1] E. Kim and M.H. W. Chan, Nature 427 225 (2004). 

[2] E. Kim and M.H. W. Chan, SCIENCE 305 1941 (2004). 

[3] G.V. Chester, Phys. Rev. A. 2,256 (1970). 

[4] T. Leggett, Phys. Rev. Lett. 25, 1543 (1970). 

[5] W.M. Saslow, Phys. Rev. Lett. 36, 1151 (1976). 

[6] A.F. Andreev and I.M. Lifshitz, Soviet Phys. JETP 29, 

1107 (1969). 

[7] N. Prokofev and B. Svistunov, [cond-mat t /0409472 

[8] R.A. Guyer, Phys. Rev. Lett. 26 174, (1971). 

[9] M.W. Meisel, Physica B 178 121, (1992). 



[10] B. A. Fraass, P. R. Granfors, and R. O. Simmons, Phys. 
Rev. B 39, 124 (1989). 

[11] B. Chaudhuri, F. Pederiva and G. V. Chester, Phys. Rev. 
B 60, 3271 (1999). 

[12] D. M. Ceperley and B. Bernu, |conTTmat /0409336 

[13] In these variational calculations shown in Figs.lc, and 
Id, for given mean fields a — b, we first solve the mean 
field Hamiltonian of (1) to obtain the lowest energy state 
| *(a)>. We then calculate E(a) = (*(o) | H \ 9(a)). 



